Hydrogenic Spin Quantum Computing in Silicon: a Digital Approach 
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We suggest an architecture for quantum computing with spin-pair encoded qubits in silicon. 
Electron-nuclear spin-pairs are controlled by a dc magnetic field and electrode-switched on and off 
hyperfine interaction. This digital processing is insensitive to tuning errors and easy to model. 
Electron shuttling between donors enables multi-qubit logic. These hydrogenic spin qubits are 
transferable to nuclear spin-pairs, which have long coherence times, and electron spin-pairs, which 
are ideally suited for measurement and initialization. The architecture is scalable to highly parallel 
operation. 

PACS numbers: 03.67.Lx, 72.25.Dc 
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A quantum computer comprising many two-level sys- 
tems, or "qubits," exhibits coherent superpositions (the 
incompatibility of certain observables) and entanglement 
(strong correlations between qubits). These quantum 
features may be harnessed to solve problems which are 
essentially impossible for a classical computer, such as 
the factorization of large integers or the simulation of 
many-body quantum systems jlj. Solid state implemen- 
tations stand to benefit from the rapid advances in semi- 
conductor electronics and are potentially scalable to large 
arrays of qubits controlled by gate electrodes. Donor nu- 
clear spins in silicon are especially good solid state qubits 
because of their long coherence times. They can in princi- 
ple be controlled by hyperfine-tuned magnetic resonance 
techniques and coupled by the electron exchange interac- 
tion when carefully tuned surface gate voltages properly 
position the donors' electrons Q. However, this "ex- 
change mediation" is restricted to nearest neighbor in- 
teractions and is extremely difficult to control [3, Q ; the 
coupling strength is very sensitive to the electrons' posi- 
tions ^exhibiting rapid oscillations due to Si band struc- 
ture 0, 0. Precise tuning of the hyperfine interaction 
will also be difficult. In this paper we present an alterna- 
tive donor spin architecture which tolerates tuning errors 
and overcomes nearest neighbor restrictions. 

Our proposal relies on the "encoding" of each logical 
qubit, a |0) + (3 |1), in the J z = subspace of a pair of 
spins: |0) = (|T4) - |4T»/V2 and |1) = (|U) + |4T»/V2. 
Encoding often results in reduced constraints on com- 
puter design 0, El ■ When the two spins are donor nuclei 
the qubit benefits from their long coherence times. On 
the other hand, measurements are facilitated when the 
two spins are electrons 0, . Following Levy, who pro- 
posed Heisenberg-only quantum computing with distinct 
magnetic moments in a static magnetic field 0, EI Il2j , 
we will show that when the two spins are an electron 
and its donor nuclear spin ("a hydrogenic spin qubit") 
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FIG. 1: Schematic of the proposed architecture. Each qubit is 
encoded in the spins of an electron and its donor nucleus. "A- 
gates" above donor sites switch the electron-donor overlap, 
and thus the hyperfine interaction, while "S-gates" shuttle 
electrons from donor to donor. "Bit trains" of voltage pulses 
control the computer. 



the qubits are easier to control and can be coupled, well 
beyond their nearest neighbors, with electron shuttling. 

In the hydrogenic spin qubit the electron and donor nu- 
clear spins are coupled by the hyperfine interaction. The 
ground state coupling for P donors in Si, Ha = Aa e - a n , 
is ideally suited to quantum computing because it is a 
quadratic, and thus insensitive, function of small ex- 
ternal electric fields (perturbations from zero field); its 
strength, A = 121.517±0.021 neV [HI, is determined by 
the electron-donor overlap, ^(O)! 2 . Here a = [a x ,<r y ,a z ) 
are the Pauli operators, labeled by the spin on which 
they operate. As depicted in Figure ^ we can use a 
surface "A-gate" voltage to draw the electron off the nu- 
cleus, effectively switching off the coupling (Ha — > 0) to 
a regime which is similarly insensitive to tuning errors. 
We therefore propose a digital approach [14j, in which 
the interaction is only on or off for each clock cycle as 
determined by "bit trains" of voltage pulses. The hyper- 
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invariant subspaces and magnetic energy levels 



FIG. 2: Entangling qubits e\n\ and e2fi2. S-gates displace 
ei and shuttle e\ to the vicinity of 712 • The A-gate above n.2 
then applies hyperfine interaction, generating a partial ei— n.2 
spin swap. 



fine control generates the electron-donor (e-n) spin swap 
|0) + |1) <-> |0) — |1) and we augment this with a globally 
applied static magnetic field, which generates |0) <-» |1). 
For 0(1 mT) fields the two generators are of comparable 
strength and an alternating series of interactions imple- 
ments single qubit logic in direct analogy with Euler's 
theorem for constructing an arbitrary rotation from a se- 
quence of rotations about distinct axes Q . 

Electron spin coherence distances of over 100 /im have 
been demonstrated yj|, so single electron shuttling 
to remote donor sites is a good candidate for enabling 
two-qubit interaction. As shown in Figure [21 arrays of 
"S-gate" electrodes between qubits are thus used to shut- 
tle individual electrons from site to site. Two qubits be- 
come entangled when the hyperfine interaction is applied 
between the electron of one qubit and the nucleus of an- 
other. This is analogous to ion-trap proposals in which 
ions, and thus their quantum information, can be trans- 
ported from one local trap to another [17lll8l |. This trans- 
port is considerably more efficient than a bucket brigade 
series of nearest neighbor interactions and can circum- 
vent misbehaved donor sites. 

The evolution of the electron and donor spins is de- 
scribed by their Hamiltonian, 
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The second term, Hb, sums the contribution from all 
donors and their electrons, with respective magnetic mo- 
ments g n pN an d 9eMS; i n the vertical magnetic field B 
assumed parallel to a (100) lattice plane. It augments the 
hyperfine contact term, Ha, which is a sum of interac- 
tions between electron-donor pairs. Interaction between 
the ith electron and the jth donor is either off (Ay = 0) 
or on (Aij = A). We assume instantaneous switching 
and neglect the hydrogenic spin-orbit and dipole-dipole 
interactions (which are zero for the ground state and for 
sufficiently large r but finite in between) as well as any 
randomness in the contact strength during the switch. 
For P donors in Si the ground and first excited orbitals 
are separated by s» 15 mcV; a more realistic adiabatic 
switch takes 0(3 ps) which is fast compared to the hy- 
perfine interaction. Any remaining hydrogenic spin-orbit 
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FIG. 3: Magnetic energy levels and invariant subspaces of 
a two-qubit computer. Flipping a single electron or single 
nuclear spin changes the energy by AE e or AE n respectively 
and takes the state to another subspace. Within an invariant 
subspace, simultaneous electron and donor spin flips change 
the energy by AE r = AE e + AE n . 



and dipole-dipole effects are coherent and can in princi- 
ple be compensated by sophisticated control sequences 
or pulse shaping 0, |2(|, although we do not consider 
them here. Similarly, we neglect the spin-orbit effect at 
the interface [2l| because, for controlled shuttling of indi- 
vidual spins in Si, it is small, coherent, and, with further 
research, characterizable and correctable. 

The state space of spins is decomposable into invari- 
ant subspaces labeled by the z component of the total 
spin; up and down spins are stationary states of Hb while 
electron-donor spin swaps, generated by Ha, preserve the 
number of up vs. down spins. Within each invariant sub- 
space flipping an electron spin, which changes the energy 
by AE e = 1Bg e [iB has a compensatory nuclear spin flip, 
which changes the energy by a further AE n = 2Bg n ixjq, 
and the magnetic energy splittings are thus integer mul- 
tiples of AE r — AE e + AE n . Transitions between sub- 
spaces require the flipping of one spin or the other and 
thus there exist nonresonant shifts AE e and AE n be- 
tween subspaces. As a specific example Figure [21 shows 
the magnetic energy levels and invariant subspaces of a 
two-qubit computer. 

It is desirable to generate pure hyperfine evolution even 
though the magnetic field is, in fact, always present. We 
make use of the Trotter formula 
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and compose a finite duration, t, of hyperfine evolution 
with a large number, a, of short At = t/a steps of hy- 
perfine and magnetic evolution corrected, on the fly, by 
time-reversed At/ 2 steps of solely magnetic interaction. 
Although magnetic and hyperfine steps do not commute, 
the remaining error of each step after correction, try 
a variant of the Campbell-Baker-Hausdorff formula yj, 
is C(At 3 ) and increases with magnetic field (the non- 
commutivity, [Ha,Hb], scales with B); we can achieve 
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good fidelity with sufficiently short At steps and a suffi- 
ciently weak field. 

Within each invariant subspace the time-reversed mag- 
netic steps are achieved by incomplete periods of mag- 
netic evolution. A full period is determined by the energy 
splitting: Tb = h/AE r (see Figure^ . We need only wait 
for Tb — At/2 to achieve the magnetic correction step. In 
analogy with magnetic resonance techniques we thus pro- 
ceed by resonant stepping; for each period of magnetic 
evolution there is a short step of Ha + Hb ■ The result 
is true hyperfine evolution up to relative phase shifts be- 
tween invariant subspaces. 

The use of digital bit trains from a pulse pattern gen- 
erator considerably simplifies the timing of these opera- 
tions. For example, we divide the fixed hyperfine period, 
T A = h/4A = 8.50847 ns, into 96 clock cycles by setting 
the frequency at / = 11.2829 GHz; given this frequency 
we then divide the magnetic period Tb into 256 clock 
cycles by choosing a field strength of B = 1.57171 mT. 
Within an invariant subspace, generating pure hyperfine 
evolution is now as simple as turning off certain A-gate 
voltages for 2 clock cycles out of every 256. 

The encoded qubits reside in the J z = invariant sub- 
space. We can thus construct logic operations from finite 
<t> pulses of magnetic evolution, (B,(ji) = e - lH -B<t>T B l~h ^ an( j 
9 pulses of pure hyperfine evolution, (A, 9) = e — lH A8T A /h ^ 
implemented with resonant hyperfine stepping. The 
Controlled-Not (CNOT) operation, which performs a log- 
ical NOT operation on a second qubit contingent on the 
state of a first, can, for example, be implemented in the 
following manner: 

CNOT = (ii ® Z 2 )N(L 1 ® Z 2 )\ 

in which single qubit operations, 

(Li ® Z 2 ) = (B, y)(4 n + A m> *)(A u , |)(B, |), 
augment an entangler, 

N = (A 12 +A 21 , ^)(B, J)(A 21 ,7r)(B, ^)(A 12 +A 21 , J). 

This construction refines Levy's original fll| but may not 
be optimal. Figure 0] depicts the actual sequence of A- 
gate voltages that implements the entangler, N. 

Shorter hyperfine steps and a weaker magnetic field re- 
duce the errors. However, commercially available pulse 
pattern generators are limited to approximately 12 GHz 
(hence our choice of / = 11.2829 GHz). Furthermore, the 
preponderance of magnetic periods (one for each small 
hyperfine step) means that a computation slows with 
weaker field. There is thus a trade-off between fidelity 
and speed. Our choice of B = 1.57171 mT yields a 
complete spin swap (the architecture's fundamental pro- 
cess) in 0.57/xs. When ideally implemented with reso- 
nant hyperfine stepping, its expected error (defined to 
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FIG. 4: Operations time line for an entangler. Between shut- 
tling operations represented schematically by an "S," the dia- 
gram depicts a sequence of on-or-off A-gate voltages and their 
duration in clock cycles. 

be the average probability of incorrectly transforming an 
initial, arbitrary, two-qubit basis of states) is less than 
2.1 x 10~ 7 . The CNOT is our most complicated gate and 
can be ideally implemented with an expected error of at 
most 0.9 x 10~ 6 in 3.22 /zs. 

It is unrealistic to presume exact values for the fre- 
quency, field, and hyperfine strength. There may also 
be variations of hyperfine and/or field strength from one 
donor site to the next. Indeed, although isotope purifica- 
tion can remove most Si 29 from the crystal, the remaining 
impurities cause field variations (although these fluctuate 
so slowly that spin-echo techniques may be applicable). 
Another complication is that the Lande factor for the 
electron, g e , could vary by as much as 10~ 3 between the 
donor and the Si-barrier interface [l3j |. 

We have studied the sensitivity to these parameters 
by the explicit simulation of canonical one- and two- 
qubit logic gates. The threshold theorem Q for quantum 
computation concludes that efficient quantum comput- 
ing, obtained with error-correction techniques, is possible 
when logic gate errors are less than 10~ 5 . We found that 
this threshold is obtainable with relative variations in 
frequency, field, and hyperfine strength as large as 10~ 5 , 
10~ 5 , and 5 x 10~ 4 , respectively. The sensitivity to lo- 
cal variations in these parameters is approximately the 
same. The fidelity is comparatively insensitive to the hy- 
perfine strength because our gate compositions are pre- 
dominantly magnetic. Finally, the architecture can tol- 
erate 5 x 10~ 3 variations in g e between the donor and the 
interface. 

A 7r pulse of hyperfine interaction, (A, 7r), between two 
qubits generates a complete spin swap between the elec- 
tron of one qubit and the donor of the other. Considered 
as a switch to a new encoding scheme, this hyperfine 
"data bus" transfers one qubit into a nuclear spin-pair 
and the other into an electron spin-pair. For example, 
an en data qubit, with the use of an e^n^ "ancilla," can 
be transferred, by resonant hyperfine stepping, into an 
riAn nuclear spin-pair qubit. Retrieval simply requires 
another 7r pulse to repeat the spin swap. 

The relatively weak nuclear magnetic moment gives 
the nuclear spin a long dccohcrcnce time which makes 
the nuclear spin-pair qubit a natural quantum memory. 
Furthermore, if the data and ancilla were unentangled be- 
fore the swap then the data (now encoded in the nuclear 



4 



e 1 n 1 e 2 n 2 





\ 


SET 


triplet 






|0)«|1) 


/ 




} 



singlet 







SET 


/ 




\ . 

singlet 


\ 



4, 

I junk> 



FIG. 5: Qubit Initialization and Sorting. A singlet outcome 
is immediately convertible into |0) while the triplet outcome 
can be recycled through a sequence of operations into another 
chance for a useful singlet. 

spin-pair) and ancilla (now encoded in the electron spin- 
pair) remain unentangled, so decoherence or collapse of 
the electron spin-pair will not degrade the memory (the 
qubit 's transfer succeeds even when the ancilla is outside 
its logical subspace; relative phases developed between 
invariant subspaces, by resonant hyperfine stepping, are 
absorbed solely into the ancilla). 

The data qubit can, alternatively, be transferred into 
an electron spin-pair to facilitate measurement by vari- 
ous proposed methods to distinguish singlets and triplets. 
For an electron spin-pair known to reside in the log- 
ical subspace, these are effectively |0) = \singlet) vs. 
|1) = [triplet, S z — 0) projective qubit measurements. 
For example, a Single Electron Transistor (SET) is capa- 
ble of very sensitive charge configuration measurements; 
above a donor it can detect electrode driven charge den- 
sity fluctuations associated with the electron spin-pair 
sing let 0. Alternatively, in a quantum dot the electrons' 
spin determines the tunneling of spin-polarized currents 

After measurement the collapsed electron spin-pair can 
be transferred back into an electron- donor pair via an- 
other spin swap. This provides a way to initialize the 
computer at high temperature (e. g., 1 K). Readout col- 
lapses an electron spin-pair into a singlet or triplet. The 
singlet outcome, |T ei |e 2 ) — |leiTe 2 )) is immediately con- 
vertible, via a spin swap, to |0). The triplet outcome, 
lteite 2 ), |T ei | e2 > + lJ-eiTe 2 ), or || ei | e2 ), can be recycled, 
as depicted in Figure [5] through a single qubit |0) «-> |1) 
operation sandwiched between spin swaps, for another 
chance to obtain a useful singlet. (This cascaded mea- 
surement prevails despite relative phases developed be- 
tween invariant subspaces.) At high temperature 50% 
of the electron-donor pairs will obtain |0), and by elec- 
tron shuttling the successful 50% can be "pooled" into 
the working part of the computer in analogy with Kane's 
original proposal for on-chip spin refrigeration |22j . 

Hydrogenic spin qubits and coherent single electron 
shuttling enable a silicon-based quantum computer fea- 
turing digital hyperfine control insensitive to tuning er- 



rors, a long-lived nuclear spin memory, a projective read- 
out scheme, and qubit refrigeration in which 50% of the 
qubits can be initialized at high temperature. The com- 
puter is scalable to highly parallel operation because dig- 
ital shuttling of electrons overcomes nearest neighbor re- 
strictions. Finally, donors can be irregularly spaced and 
far apart, allowing for large gate electrodes, and malfunc- 
tioning donor sites can be diagnosed and avoided. These 
many benefits motivate further research on the coherent 
shuttling and measurement of electron spins, extremely 
pure Si fabrication, encoding and error-correction tech- 
niques, optimal control sequences, and the spin-orbit 
and dipole-dipole interactions during realistic electrode 
driven switching and shuttling. 

We are grateful for helpful discussions with S. 
Lomonaco. 
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